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We propose a selective dynamical decoupling scheme on a chain of permanently coupled qubits 
with XX type interactions, which is capable of dynamically suppressing any coupling in the chain 
by applying sequences of local pulses to the individual qubits. We demonstrate that high-fidelity 
single- and two-qubit gates can be achieved by this procedure and that sequences of gates can be 
implemented by this pulse control alone. We discuss the applicability and physical limitations of our 
model specifically for strongly coupled superconducting flux qubits. Since dynamically modifying 
the couplings between flux qubits is challenging, they are a natural candidate for our approach. 
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I. INTRODUCTION 

Current implementations of qubits are typically ei¬ 
ther well isolated from noise, but difficult to couple, or 
strongly coupled, but difficult to isolate. To achieve the 
goal of building a working quantum computer for the 
tasks of quantum simulation and, eventually, quantum 
computation, we require an architecture that is capa¬ 
ble of strongly coupling qubits to implement fast multi¬ 
qubit gates, but that can also isolate qubits from each 
other and the environment when no gate operation is 
performed. The trade-off between strong coupling and 
isolation should be optimized to maximize the ratio be¬ 
tween decoherence time and gate operation time. 

In quantum optics, extensive work has been done us¬ 
ing trapped ions or atoms as qubits, and scalable archi¬ 
tectures that can trap and address a large number of 
qubits simultaneously exist [I|. These qubits feature ex¬ 
cellent coherence times, yet the implementation of two- 
qubit gates in these architectures is still a topic of ongo¬ 
ing research, although recently promising proposals were 
made in this regard M- 

Likewise, we have seen s^nificant progress in solid 
state qubit architectures (5l-l7|, and there exist promis¬ 
ing candidates for scalable systems. Gate-defined spin 
qubits d, d feature excellent coherence properties [I^ , 
but coupling two qubits remains a challenge despite pro¬ 
posals for efficient coupling [n|, [T^ ■ For superconduct¬ 
ing qubits, both indirect coupling via a resonator [I^ 
and direct capacitive coupling of detuned qubits [l^ have 
been demonstrated and are comparatively easy to real¬ 
ize. While in principle superconducting qubits can be 
coupled very strongly, the small anharmonicity of some 
designs limits the possible coupling strength [l5j. Very 
recently dynamically protected superconduc ting qubits 
have been proposed [1^ and demonstrated [l3l, which 
rely on a parametrically driven cavity [l8l | in the quan- 
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turn regime (l9l-l^. In this design the qubits are well 
protected, but direct coupling of two qubits might re¬ 
main a challenge. 

Good coherence properties were achieved for flux 
qubits [1^, a particular type of superconducting qubit 
with a very large anharmonicity. This anharmonicity al¬ 
lows them to be strongly coupled [2^ , which makes them 
particularly interesting for the implementation of fast 
two-qubit gates. Larger circuits containing many we akly 
coupled flux qubits have already been demonstrated . 
However, their tunability is limited by the need for an op¬ 
timal operating point, which makes it difficult to isolate 
the qubits when no gate operation should be performed. 

In this article, we attempt to overcome the isolation 
problem of flux qubits and other similarly strongly cou¬ 
pled qubit systems by an alternative ansatz. We study 
a qubit chain with always present nearest-neighbor cou¬ 
plings and make use of a pulse generator to exert ex¬ 
ternal control on the qubits with the aim of suppressing 
unwanted qubit couplings. We demonstrate in numerical 
simulations that this simple pulse control enables us to 
implement a sequence of entangling gate operations on 
the qubit chain to entangle all the qubits in the chain in 
a GHZ state [1^ with high fidelity. We thus show that 
a system of strongly coupled flux qubits may be used 
for universal quantum computation purposes without the 
need to control the qubit couplings. 

Our pulse control is based on dynamical decoupling 
, which is a generalization of techniques developed in 
the nuclear magnetic resonance (NMR) community 113- 
1^ . It makes use of external control pulses being applied 
in rapid succession to the system in question. With a 
carefully designed control sequence it is possible to elim¬ 
inate (parts of) a Hamiltonian interaction up to a cer¬ 
tain order. Dynamical decoupling has been successfully 
implemented in numerous experiments to protect qubit 
states from the effects of decoherence [3ll - [^ . For our 
purposes, we are interested in selectively decoupling only 
certain interactions between qubits while keeping oth¬ 
ers alive, a possibility proposed already by Viola et al. 
(^ . A particularly simple to handle subset of decou¬ 
pling schemes applicable to networks of qubits employs 
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only Pauli pulses to individual qubits. Several different 
construction methods for such Pauli operator schemes 
exist (37l - l4lj |. In dynamical decoupling, it is typically as¬ 
sumed as a first approximation that the applied control 
pulses are instantaneous and unitary. In our numerical 
calculations, we go beyond this approximation by simu¬ 
lating realistic pulse lengths. To deal with such bounded 
controls, advanced decoupling techniques in the form of 
Eulerian decoupling and dynamically corrected gates 
( 4 ^ 1^ have been developed, which we will make use of. 

The paper is organized as follows: In section |ll] we 
present the physical model of our qubit chain and the 
type of control we have over the system. Section |TTT] ex¬ 
plains how a two-qubit iSWAP gate can be implemented 
with the help of dynamical decoupling. Decoupling basics 
are explained and a decoupling sequence for this partic¬ 
ular task is developed. Numerical results for the achiev¬ 
able fidelity are presented. In section HVl we then look at 
how to implement single-qubit gates with high fidelity, 
where numerical simulations were conducted to verify the 
achievable fidelities. Finally, in section we introduce 
the CNS gate and use this gate to entangle all the qubits 
in our chain in a GHZ state. We calculate numerical re¬ 
sults for the achievable GHZ state fidelity for different 
numbers N of qubits and also look at how much of an 
impact disorder has on the fidelity. 


II. THE COUPLED QUBIT SYSTEM MODEL 

We consider a system of N qubits in a chain with 
nearest-neighbor couplings described by the Hamiltonian 

N N-1 

where the are the Pauli operators applied to the i-th 
qubit, and are the qubits’ eigenenergies. The coupling 
between the qubits is assumed to be uniform and charac¬ 
terized by the coupling strength g. This model is strongly 
inspired by a system of coupled flux qubits [1^, however, 
alternative qubit designs exist which are also described 
by this Hamiltonian. Additionally, in our model there is 
a pulse generator with frequency oj which can exert ex¬ 
ternal control on the qubits, and in the case of flux qubits 
is implemented as a microwave emitter. It is described 
by the control Hamiltonian 


frame and with the rotating wave approximation, the sys¬ 
tem and control Hamiltonians equal 


N-l 
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. N 

Hr{t) (cos{ip,(t))a["^ -tsin((/j,(t))CT^*^) . 


(3) 


The Ai = ei — OJ indicate the detuning between the in¬ 
dividual qubits’ eigenenergies and the frequency of the 
driving field and should ideally be zero for our purposes. 
If the eigenenergies are different, then we have disorder, 
which can disrupt the gate operations we intend to im¬ 
plement in the following. However, as we will see, our 
approach is robust to disorder due to our use of decou¬ 
pling, as long as the do not become too large. 


III. IMPLEMENTING THE TWO-QUBIT 
iSWAP GATE BY SELECTIVE DECOUPLING 


The coupling between the qubits according to ([3]) is 
of XX type. Schuch and Siewert [i^ studied natural 
gate operations resulting from such an interaction. They 
showed that, after an interaction time T = Tr/{2g), this 
type of coupling between two qubits produces a unitary 
iSWAP gate: 


GiswAP :=exp 







-I- fJo cr. 


(4) 


This gate, like the better known SWAP gate, exchanges 
the state of two qubits, but introduces an additional 
phase on the swapped qubit states. However, in our 
model we have additional couplings to the qubits (f — 1) 
and {i + 2) as well as the disorder terms A^crg*^ and 

Ai-i-KTg*’''^^ In order to succesfully use the natural cou¬ 
plings to implement the iSWAP gate, we need to isolate 
the two qubits involved in the gate operation. Tradi¬ 
tionally, we would thus require switching off any interac¬ 
tions which are not currently needed, but this process is 
complicated and often limits the achievable interaction 
strength g. Instead, we will employ dynamical decou¬ 
pling to suppress the effects of individual couplings as 
needed. 


A. Dynamical decoupling basics 


N 

Hco{t) = ^ /i(t)cri*^ cos(wt -f (2) 

i=l 

and is governed by the pulse amplitudes fi{t) and phases 
which can be controlled for each qubit individually. 
It is convenient to switch to a rotating frame by trans¬ 
forming to the interaction picture given by the unitary 
operator Uai{t) = exp(ia; cr3*^t/2). In the rotating 


In dynamical decoupling, the natural evolution of the 
V-qubit chain under the acting Hamiltonian H is modi¬ 
fied in a controlled fashion by the external control Hamil¬ 
tonian Hc{t). In our case, the pulse generator will be 
activated periodically at times tj for a short time tp to 
implement a sequence of pulses 


Pj = +tp) =T exp 


ptj ~\~tp 


dtHcit) 


(5) 
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where T denotes the Dyson time-ordering operator and 
Pj is a unitary operator representing the j-th pulse of 
the sequence. However, the implementation of the pulse 
is disturbed by the acting Hamiltonian H, so that we get 
an imperfect pulse of the form 


Pj = 'P exp y- 


dt{Hcit) + H)^ = Pjp]pj 


ptj -\-tp 


dtUl{tj, t)HUc(tj , t) 


( 6 ) 


If we assume that after each decoupling pulse there is a 
time T of free evolution under the Hamiltonian H, then 
by introducing the unitary operators pj = pjpj-i ■ ■ -po, 
the resulting time evolution U{t) after M pulses can be 
written as 



U{Mt) =pMe-^^^pM-i ■ ■ 

=9M pM-ij 

. . . g-Do^iSOg-iSoffso-r 

It is customary to enforce the cyclic condition pM = Po = 
1 by an appropriate choice of the decoupling pulses pj. 
We can now define an average Hamiltonian H which leads 
to the same time evolution after the time Mr, i.e., 


FIG. 1: A Eulerian path decoupling sequence for a single 
qnbit. I,X,Y,Z correspond to the Pauli operators 1, cri, CT 2 
and as. The vertices represent the decoupling operators pj, 
the directed edges denote the transitions between the pj due 
to the decoupling pulses pj. 

B. Decoupling an individual qubit 

As a first step, we will discuss decoupling of a single 
qubit on the chain with the goal of freezing the evolution 
of that qubit’s state. For a single qubit, there exists a 
particular decoupling scheme 


U{Mt) = 

By performing a Magnus expansion the average 

Hamiltonian H is expanded in powers of the pulse dis¬ 
tance T, i.e., 


H = H 


[ 0 ] 


-b H 


[ 2 ] 


(9) 


where the lowest order is found to be 

M-l 


H 


[ 0 ] 


3=0 


H +-^ 
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( 10 ) 


Here, is the lowest order of the Magnus expansion of 
the error operator which is given by 

dtUl{t)HUj,^{t). ( 11 ) 

do 

Our goal is to selectively remove couplings between 

specific qubit pairs in the lowest order of the aver¬ 
age Hamiltonian and to keep all others, while simultane¬ 
ously suppressing the effects of the disorder terms A^o-g*^ 
and the pulse errors We call a set of M operators 
{9j}fLo^ a decoupling scheme if it fulfils this purpose. 
Note that the higher orders of H are typically non-zero 
and remain as errors. 


{50 = 1, pi = cn, P 2 = CT3, 53 = CT 2 } , (12) 

which has the property that for any traceless Hermitian 
operator P acting on the subspace of the qubit, 

3 

Y.plPp^=0■ (13) 

i=0 

It can be implemented solely with the help of cti and a 2 
pulses to the qubit (meaning tt pulses around the X or 
Y axis, respectively), by the sequence 

Pl=yi 0-2, CTl, ‘’’2, Tl n A\ 

Po =1 ->• cn - 0-3 -)> (72 -1- (14) 

If we insert this decoupling scheme into (flUl) , it will elim¬ 
inate all parts of the Hamiltonian H acting on the qubit 
in the lowest order, effectively decoupling the qubit from 
the rest of the chain. Unfortunately, the pulse errors 
depend on the particular pulses pj and are thus not 
eliminated by this decoupling scheme. 

There is a trick to construct a decoupling sequence 
from the scheme in Eq. m which not only eliminates 
H, but also the pulse errors in (1101) . This method 
is called Eulerian path decoupling and proceeds as 
follows. A graph is constructed from the scheme m 
where the scheme operators pj are taken as the vertices 
of that graph and a directed edge is placed between two 










4 


operators Qa, gb if cFiga = Qb or (723a = 3b, up to a phase 
factor. A Eulerian path through this graph is a path 
which visits every edge of the graph exactly once. A 
particular Eulerian path is depicted in figure [TJ which 
results in the following decoupling sequence: 

(72 <yi (72 

go =1 —s> CTi —> ao —)> 0-2 — 

1 ->■ tT2 -CT3 -CTl - > 1. (15) 

It corresponds to two consecutive applications of the orig¬ 
inal scheme (ini), but with different orders of the scheme 
operators gj . If inserted into (ITOl) , we get 

E d +)(«;"'+4"')) »=o. 

{l.CTi,cr2,cr3} 

_( 16 ) 

The remaining orders of the average Hamiltonian H are 
of order 0{\\H\\^t) + (!l(||$j|p). 


C. Selective decoupling on the qubit chain 

The decoupling sequence discussed in section ITlI B I can 
isolate a single qubit from the chain. We need to ex¬ 
tend this sequence to the whole chain in such a way that 
we can selectively decouple only certain qubit couplings 
while keeping others alive. In [4l| we found that if we 
extend the original decoupling scheme (fl^ to two qubits 
in the following way, 



FIG. 2: The pulse sequence used to implement the iSWAP 
gate. This figure shows the pulse sequence used for both of 
the gate qubits, where blue signifies a pulse in X direction 
and red signifies a pulse in Y direction. Neighbouring qubits 
use the same pulse sequence, but with X and Y swapped if 
their coupling is to be eliminated. 


qubit, we would use the following pulse sequence: 


_ ^(1)„(2)„(3)^(4) (5) 


Pi = 


{T-i ^ ^ 


(l)^(2)^(3) (4) (5) 


P2 — ^2 ^2 
P3 = PO 
P4=Pi¬ 


rn 


While the previously explained extension allows us 
to selectively decouple certain qubit couplings from the 
Hamiltonian H in m, it does not eliminate the pulse 
errors . A simple modification to (IT^ sees us adding 
additional pulses 


5o = 1 0 1, gi = (71® (Ti, 

32 = 0-2 0 0 - 2 , 33 = 0-3 0 ( 73 , ( 17 ) 

then in OT it will keep the Heisenberg-type coupling 
terms between these two qubits intact while still elimi¬ 
nating the disorder terms Aia^\ On the other hand, if 
instead we choose 

30 = 1 0 1, gi= <7l® ( 72 , 

32 = (72 0 ( 71 , 33 = (73 0 ( 73 , (18) 

then the Heisenberg-type couplings in m are elimi¬ 
nated between these two qubits. The first scheme can 
be realised by applying the pulse sequence XYXY on 
both qubits simultaneously, while the latter one applies 
XYXY to one qubit and YXYX to the other. 

It is straight-forward to extend these schemes to the 
whole qubit chain. On each qubit, we apply alternat¬ 
ing XY pulses. Neighbouring qubit pairs whose inter¬ 
action should be kept alive will employ XY pulses in 
the same order, whereas between qubit pairs whose in¬ 
teraction should be decoupled, we use alternating pulse 
sequences. For example, if we wanted to protect the in¬ 
teraction between the first and last qubit pairs on a 5- 
qubit chain, but eliminate the couplings with the middle 


P5=P4, P6=P3, P7=P2, P8=Pl, (20) 

with the effect that on each individual qubit we now have 
a Eulerian path decoupling sequence as in (fTSl) . without 
changing the effects of the sequence on H in the lowest 
order OT. This is the final decoupling sequence which 
we will use to implement our quantum gates, and we will 
see in numerical simulations that it produces sufficiently 
high fidelities. However, we should point out that, due 
to the extension of the sequence to the whole chain, the 
pulse errors are not fully eliminated in the lowest order 
even with the Eulerian path modification. The reason 
is that the original schemes (CZD and (1151) only eliminate 
certain Hermitian operators on the two-qubit subspace, 
but not all of them. A more sophisticated approach is 
outlined in which eliminates errors completely (in 
the lowest order), however, it requires 64 pulses instead 
of 8 and thus has a significantly longer implementation 
time. 


D. The iSWAP gate and physical limitations 

We have all the necessary prerequisites to implement 
iSWAP gates on our qubit chain. The procedure is sim¬ 
ple: over the implementation time T = tt/( 2g), we apply 
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the sequence of eight pulses developed in section IIII Cl 
Due to the selective decoupling, we can implement sev¬ 
eral iSWAP gates in parallel, provided that any two gates 
do not share a gate qubit. Both cti and U 2 pulses can be 
implemented with our pulse generator as tt pulses around 
the X 01 Y axis. For cti pulses, the phase is chosen 
to be 0, for (J 2 pulses it is chosen as 7r/2. The amplitude 
fi{t) can be any smooth function with the condition 



The pulse implementation time tp should be made as 
small as possible to reduce the pulse errors. However, 
there are some fundamental obstacles which prevent us 
from making tp infinitely short. For one, a physical pulse 
generator will have limitations on how quickly it can steer 
the pulse amplitude and on the maximal achievable pulse 
amplitude, which in turn limits the minimal pulse dura¬ 
tion. Additionally, the rotating frame Hamiltonian in 
Eq. ([3]) was derived in the rotating wave approximation. 
In order to ensure validity of this approximation, we re¬ 
quire 1 <C 2ujtp. Another fundamental problem is the 
fact that many physical implementations of qubits are 
only approximately two-level systems. If we probe the 
physical system hard enough, which in our case means if 
we choose tp —>■ 0, eventually we will excite higher states 
or invoke additional interactions and thus invalidate our 
two-level approximation. 

With that in mind, let us look at what kind of pulse 
duration we would have to achieve to actually imple¬ 
ment the iSWAP gate with high fidelity. Given the im¬ 
plementation time T = nli^g) of the iSWAP gate and 
the necessity to implement a series of eight pulses during 
that time, the upper limit for the pulse time is given as 
tp < ’’■/(Ibg). In our simulation, we used pulse times 

tp € [7r/(16g),7r/(32g),7r/(48g),7r(645),7r/(965)]. 

We simulated a qubit chain of varying length with = 0 
and implemented the iSWAP gate in the middle of the 
chain. We used Gaussian pulse shapes for the decou¬ 
pling pulses, and Fig. [2] depicts the pulse sequence used. 
We simulated the time-dependent Schrodinger equation 
for the full pulse sequence and calculated the emerging 
state of the qubit chain, where we then traced out all 
of the qubits except for the two gate qubits. The re¬ 
sulting state p was then compared to the expected state 
1^') = GiswApl^'in) by means of the state fidelity [4^ 

F(T) = |(4/|p|vI/)|. (22) 

As initial states I'l'in) we used all four basis states |00), 
|0I), |I0) and |II) and took the average over the achieved 
fidelities. The remaining qubits were always prepared in 
the state |0). The average fidelities depending on the 
pulse duration tp are given in table HI The results were 
virtually independent of the number of total qubits N in 
the chain. We can see that even for the longest possible 


tp = 7r/(16p) 

V(32g) 

tt/ (48^) 

tt/(645) 

■K/{96g) 

0.9922 

0.9979 

0.9990 

0.9994 

0.9997 


TABLE I: Numerical simulation results for the achievable fi¬ 
delity of the ISWAP gate, depending on the pulse duration 

tp . 


pulse duration tp = 7r/(16g), the gate fidelity is quite 
good. 

Are these pulse durations realistic? Let us consider 
as a concrete example two superconducting flux qubits. 
Flux qubits with always-on couplings of the order of 
g ^ 500 MHz were realized in , which would allow 
for a fast implementation of the iSWAP gate. Addi¬ 
tionally, flux qubits feature a rather large anharmonicity, 
meaning that the higher energy levels after the two qubit 
states are separated by a significant gap. With a typical 
splitting of about 5 GHz between the first two levels, we 
could in theory have a pulse amplitude of several GHz 
before we risk exciting the higher states. Let us assume 
that we could safely employ a maximum pulse amplitude 
/max = lOGHz. Then the achievable minimal pulse du¬ 
ration for that amplitude depends on the specific pulse 
shape. For a Gaussian pulse like we used in our simula¬ 
tions we find that for tp = 7r/(16g), the required maxi¬ 
mal pulse amplitude is /max ^ 45^. However, with the 
assumed values of g and /max for the flux qubits, we only 
achieve a ratio of /max/ff ~ 20. As a consequence, we 
would have to reduce the coupling constant by a factor 
of about 2. Alternatively, one could also look at differ¬ 
ent pulse shapes. For example, a sine-shaped pulse would 
only require /max/<? ~ 25, which is much closer. However, 
we also found in our simulations that the sine pulse per¬ 
forms slightly worse in terms of achievable gate fidelity. 
As such, there is a compromise to be made between min¬ 
imizing the gate duration T oc 1/g and maximizing the 
gate fidelity. 

Let us assume that we choose to engineer a coupling 
strength oi g = 100 MHz, which gives us some additional 
reserves and allows us to aim for a pulse duration of 
tp = Tr/{32g) « 1 ns without exciting higher states. With 
the driving field frequency cj tuned to the approximate 
qubit level splitting of 5 GHz, this pulse time is then one 
order of magnitude larger than l/(2w), so that the rotat¬ 
ing wave approximation is still valid. The implementa¬ 
tion time of the iSWAP gate is T ~ 16 ns, during which 8 
pulses need to be applied, resulting in a pulse frequency 
of 500 MHz. The requirements for our pulse generator are 
ambitious, but not impossible. Even more encouragingly, 
in recent experiments flux qubits have been demonstrated 
with decoherence times of the order of 10 ps [23,I3. This 
means that the gate operation time is almost three orders 
of magnitude faster than the decoherence time, making 
this procedure viable for flux qubits. Other implementa¬ 
tions of the basic model from Sec. El may impose very 
different limitations. 

In the interest of maximizing the fidelity, we should 
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also point out that there exist more sophisticated pulse 
shapes than Gaussian or sine-shaped pulses. Some of 
these pulse shapes were specifically engineered to reduce 
their own error (see, e.g., for a review of NMR pulse 
shapes or for a more recent design) jOT are less likely 
to excite higher states in the system Both of 

these properties might help to improve the gate fidelity 
further. However, specifically with the self-correcting 
pulse shapes, the price to pay is typically a significantly 
higher ratio /max/5 to implement a particular pulse in the 
same time span. Thus the qubit interaction strength g 
would have to be reduced even further, meaning that de¬ 
coherence becomes a potentially larger concern. Which 
pulse shape is the most adequate depends on the spe¬ 
cific needs of a particular experiment. In our numerical 
simulations, Gaussian shaped pulses proved to provide a 
suitable compromise between achievable fidelity and re¬ 
quired maximal pulse amplitude. 


IV. IMPLEMENTING HIGH-FIDELITY 
SINGLE-QUBIT GATES 

In addition to the two-qubit iSWAP gate, we will also 
need to be able to perform single-qubit gates on the indi¬ 
vidual qubits. For the implementation of the single-qubit 
gates, we will again make use of the pulse generator. This 
means that the available gate operations are given by the 
unitary propagator of Hc(t) in ([3|). In particular, we can 
implement rotation operations around the X and Y axes, 

Ry{(j)) = (23) 

which can be realised by choosing the phase ipt appropri¬ 
ately and engineering the pulse amplitude function such 
that dtf{t) = (j). However, as with the decoupling 
pulses in section llil A[ the gate operation is disturbed by 
the system Hamiltonian H, which limits the achievable 
gate fidelity. For a single-qubit gate, we typically want to 
achieve fidelities well above 0.99, which is a requirement 
to add quantum error correction later. 


A. Dynamically corrected gates with Eulerian path 
decoupling 

Fortunately, there is a way to embed a gate operation 
Q into a Eulerian path decoupling sequence such that it 
decouples the error of the gate. This technique is called 
dynamically corrected gates (DCG) and was introduced 
in [ 4 ^ . The idea is deceptively simple. Remember in 
the original Eulerian path construction in Fig. [U each 
decoupling pulse formed an outgoing edge from every ver¬ 
tex, ensuring that its error would be decoupled to lowest 
order. We can add the identity operation 1 as another 
’’generator” to this picture, which can be represented as 
loops which go out from each vertex and point back to 



FIG. 3: A Eulerian path for a dynamically corrected gate 
operation Q. 


that same node. Let us now consider that our identity op¬ 
erations are not perfect, but in fact given by / = 
carrying an error 4>/ like the other decoupling pulses. 
Then this design ensures that the error is decoupled to 
first order. Finally, let us replace the final identity oper¬ 
ation with the actual gate Q that we want to implement, 
and let us assume that Q has the same error as the faulty 
identity operations, = <!)/. The updated graph for the 
resulting decoupling sequence is depicted in figured The 
net operation of this sequence without any errors would 
be the gate Q, as intended. Furthermore, the errors of all 
occurring operations are corrected to first order by the 
Eulerian path design. 

This design hinges on the question whether we can find 
a faulty identity operation which has the same error as 
the gate Q. It was shown in that this is possible at 
least to first order of the error. Gonsider an arbitrary gate 
Q with its time propagator given by Uqii:,) during the 
implementation time tq . We can introduce a scaled gate 
Qi/2 with time propagator Uq^^^{t) = Uq{t/2), which 
obviously needs an implementation time of 2tq to im¬ 
plement the original gate Q. It can be shown that this 
scaled gate implementation carries the same error to low¬ 
est order as the faulty identity gate / = Q'I'Q with time 
propagator 


\Uq{2tg-t), tq<t<2tq. 

In our control scheme, for any of the possible rotation 
gates Ra{'~p), the gates / and Q1/2 can be implemented 
in a straight-forward manner by modifying the phase am¬ 
plitude functions For the faulty identity gate I we 
need 


nit) 


Mt), 0<t<tq, 

-M2tq-t), tq<t<2tq, 


(25) 
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FIG. 4: The pulse sequence applied to a single qubit to im¬ 
plement a dynamically corrected Rxiir/2) gate with Gaussian 
pulse shapes. Blue indicates that the pulse generator is acting 
along the X axis, red indicates a pulse along the Y axis. 


meaning that we add the negative reverse of the original 
pulse shape. For the gate Q 1/2 we need to scale both the 
time and the amplitude by 1/2, meaning 

n{t) = \h{t/2). (26) 

If our minimal gate time is given by tp, then each of the 
faulty I operations and the final gate Q will take 2tp 
to implement. As a consequence, the total duration to 
implement a single-qubit gate is 16tp. For the case of flux 
qubits as discussed in Sec. IIIIDi the operation times for 
a single gate and the iSWAP gate are comparable. 

B. Implementation and numerical simulations 

Figure 0] shows the concrete pulse sequence we are em¬ 
ploying in our numerical simulations to implement a dy¬ 
namically corrected Rxij/2) gate with the decoupling se¬ 
quence from figurelSl The blue parts indicate pulses along 
the X axis, red parts indicate pulses along the Y axis. All 
qubits in the chain are subjected to the same sequence, 
except that neighbouring qubits will have the X and Y 
pulses interchanged such that the couplings between the 
qubits are decoupled. Qubits on which no gate is im¬ 
plemented will leave the pulse amplitude set to 0 during 
the / and Q phases in the sequence. Note that several 
single-qubit gates can, in principle, be applied in parallel 
to different qubits, however not on neighbouring qubits. 
The reason is that on neighbouring qubits, the error as¬ 
sociated with the gate Q = Qi®Q 2 contains terms which 
cannot be decoupled by our decoupling scheme, and as 
a consequence the fidelity reduces significantly. There¬ 
fore, single-qubit gates on neighbouring qubits should be 
performed sequentially. 

As with the iSWAP gate, we simulated the pulse se¬ 
quence from Fig. 0] on the middle qubit of a chain with 
N qubits by simulating the time-dependent Schrddinger 
equation, then tracing out all qubits but the gate qubit. 
The remaining traced state p was compared to the ex¬ 
pected state. As input states, we simulated both |0) and 


tp = 71/(165) 

tt /(245) 

V(325) 

tt/( 4 O 5 ) 

tt /(485) 

0.99929 

0.99986 

0.99996 

0.99998 

0.99999 


TABLE II: Numerical simulation results for the achievable 
fidelity of the Rx{i^/2) gate, depending on the pulse duration 



FIG. 5: The quantum circuit to implement a CNS gate with 
the help of the iSWAP gate and a number of single-qubit 
rotations. 

|1) and took the average of the resulting fidelities. The 
results for the implementation of the Rx{tt/2) gate can be 
found in table I//I Results for different single-qubit gates 
are very similar. We can see that even for tp = Tv/iXQg) 
the fidelity is excellent. 


V. ENTANGLING THE CHAIN QUBITS WITH 
THE HELP OF A CNS GATE SEQUENCE 

In the following, we investigate how to implement an 
entangling two-qubit gate in our model. An entangling 
gate is a necessity for universal quantum computing, and 
the previously implemented iSWAP gate on its own is 
not sufficient. However, the iSWAP gate can be com¬ 
bined with a sequence of single-qubit gates to perform 
the so called CNS gate [1^ , which is a combination of a 
standard CNOT followed by a SWAP operation. The 
gate sequence depicted in figure [5] implements a CNS 
gate with the upper qubit being the control qubit. If 
the control is in the state 1, then the state of the second 
qubit is flipped. Afterwards, the states of both qubits 
are swapped. This gate is able to generate entanglement 
between two qubits. 

In [15 , Hadamard gates and rotations around the Z 
axis were used. We rearranged the gate sequence to use 
rotations around the X and Y axes instead, as these are 
the operations accessible in our model with the help of 
the pulse generator. 

We already have all the pieces of the puzzle to imple¬ 
ment the CNS gate. Given that the single-qubit gates 
must be performed sequentially due to being on neigh¬ 
bouring qubits, the CNS gate will take time + 

\12tp to implement. For our flux qubit example with 
a coupling strength oi g = 100MHz and tp = 7r/(32(;), 
this yields a time of approximately 126 ns, which is still 
a factor of 80 below the decoherence time. 

As a final experiment in this chapter, we will perform 
a sequence of CNS gates to entangle all the qubits in the 
chain. 
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N 

tp = 77 /( 165 ) 

77 /( 325 ) 

77 /( 485 ) 

77 /( 645 ) 

77 /( 965 ) 

3 

0.964 

0.989 

0.995 

0.997 

0.999 

4 

0.933 

0.982 

0.992 

0.995 

0.998 

5 

0.882 

0.974 

0.988 

0.993 

0.997 

6 

0.835 

0.967 

0.986 

0.992 

0.996 

7 

0.821 

0.962 

0.983 

0.990 

0.996 

8 

0.784 

0.956 

0.981 

0.989 

0.995 

9 

0.710 

0.947 

0.977 

0.987 

0.994 


FIG. 6: A quantum circuit to entangle all qubits in a quantum 
register in a GHZ state. In this figure, the CNS gates are 
represented by a directed CNOT gate followed by a SWAP 
gate. 


TABLE III: Numerical simulation results for the achievable 
fidelity of the GHZ state, depending on the number N of 
qubits and the pulse duration tp. 


A. An entangling sequence of CNS gates 


If we perform a CNS gate on two qubits, of which 
the first (control) is prepared in the superposition (| 0 ) + 
|l ))/-\/2 and the second in the state | 0 ), then the result¬ 
ing state is (|00) -I- |ll))/-\/2, which is an entangled Bell 
state. If we now take a third qubit, initially also in the 
state |0), and perform a CNS gate an qubits 2 and 3, then 
we get a three-qubit entangled state (|000) -I- |lll))/-\/2. 
With each additional execution of a CNS gate, we can 
bring an additional qubit into the entangled state. This 
typ e of multi-qubit entangled state is called a GHZ state 


|GHZ) 


V2 


(27) 


Let us assume that all qubits on the chain are initially 
prepared in the state |0). Then we bring a qubit in the 
middle of the chain into the superposition (|0) -I- \ l))/y/2. 
This is done by applying a Hadamard gate to it, which 
in our model we can express as an X gate followed by a 
rotation Ry (— 7 r/ 2 ). From there on we apply CNS gates 
to entangle this qubit with all the other qubits in the 
chain, where we can in fact apply CNS gates in parallel. 
A gate sequence for a 6 -qubit chain is depicted in figure 

m 

We conducted numerical simulations for this gate se¬ 
quence by calculating the resulting state Idt) by simulat¬ 
ing the time-dependent Schrodinger equation, where we 
assume that all qubits are initially in the state |0). We 
calculated the fidelity Fqhz of the GHZ state depending 
on the pulse duration tp for Gaussian pulse shapes. 


Fghz = |(GHZ|«-)|. (28) 


We simulated qubit chains of up to 9 qubits. The re¬ 
sults are shown in table IHII Given pulses which are suffi¬ 
ciently quick compared to the coupling strength g, a high 
fidelity of 0.99 for the entangled state can theoretically 
be achieved even for N = 9 qubits. However, at least in 
the flux qubit case, this would require to reduce the cou¬ 
pling strength g to the point that the full gate sequence 


will approach the flux qubit decoherence time. For the 
more realistic pulse duration tp = 7r/(32g) the achieved 
fidelities are not as spectacular, but still promising. 

It is clear that with increasing N, the fidelities will 
steadily drop. This is a consequence of the increased 
number of imperfect gate operations. Additionally, the 
longer the gate sequence, the closer we get to the de¬ 
coherence time, at which point everything breaks down. 
In order to achieve scalability, the addition of quantum 
error correction is therefore necessary. We believe that 
the demonstrated gate fidelities for single-qubit gates 
and the iSWAP gate are sufficiently high that error cor¬ 
rection is feasible. For a possible implementation, we 
would propose to extend the qubit chain model to a two- 
dimensional grid, on which we could then employ a sur¬ 
face code. The extension to the grid requires modifica¬ 
tions to the decoupling scheme, which are not trivial, 
but should be possible. Such a scenario has been accom¬ 
plished recently for Ising-type qubit couplings by De and 
Pryadko in [llltl. 


B. Influence of disorder 


The results in table m were achieved under the as¬ 
sumption that the qubits’ eigenergies are in resonance, 
meaning that the in Eq. m are all zero. Non-zero 
Ai have a detrimental effect on the achievable fidelity. 
However, our decoupling scheme offers limited robust¬ 
ness against these effects. We ran additional simulations 
where we sampled the randomly from a Gaussian dis¬ 
tribution with mean value /r = 0 and standard deviation 
a. Results of the achievable fidelity depending on cr, av¬ 
eraged over 100 runs, are plotted in Fig. [3 for a chain 
of four qubits. We can see that the drop in the aver¬ 
aged fidelity is noticeable for tp = 71 /( 165 ), but with 
faster pulses becomes negligible, at least up to the sim¬ 
ulated maximal value of ajg = 1. In a recent experi¬ 
ment with 20 flux qubits 2^ , deviations of up to 1 GHz 
were observed in the eigenenergies, which may be two 
to ten times larger than the coupling 5 , depending on 
how strongly the qubits are engineered to interact. As 
such, current experimental deviations may be larger than 
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FIG. 7: Averaged fidelity for a GHZ state achievable on a 4- 
qubit chain for different values of the pulse duration tp (Gaus¬ 
sian pulse shapes were used), when the qubit eigenenergies 
differ from each other. The Ai are randomly sampled from a 
Gaussian distribution with standard deviation a. The plotted 
results were averaged over 100 runs. 

our decoupling scheme can handle. However, we expect 
that with improved manufacturing processes the qubit 
eigenenergy discrepancies will become sufficiently small 
in the future so that the detrimental influence of the dis¬ 
order is negligible with sufficiently fast pulses. 

VI. CONCLUSIONS 

We presented a coupled qubit system modelled after 
superconducting flux qubits which is fully controlled by 
a pulse generator. The qubits are strongly coupled to 
their neighbours, and the coupling is always present. We 
demonstrated how the pulse generator can be used to 
implement both single-qubit rotations and the two-qubit 
iSWAP gate. For the implementation of the two-qubit 
gate we exploit the coupling between the qubits and use 


a Eulerian decoupling scheme to decouple the gate qubits 
from the remaining qubits in the system. The decoupling 
scheme is flexible so that several two-qubit gates can be 
implemented in parallel. The single-qubit rotations are 
realised with the help of dynamically corrected gate op¬ 
erations, which embed the gate operation into a Eulerian 
decoupling sequence. 

The efficiency of our control scheme was analyzed in 
numerical simulations, where we first looked at single 
gate applications and achieved high fidelities for both the 
iSWAP gate and the single-qubit rotations. Then a se¬ 
quence of CNS gates was simulated to entangle all the 
qubits in the chain in a GHZ state. In order to entangle 
N qubits in a GHZ state, N — 1 CNS gates are required. 
Without error correction, the GHZ state fidelity directly 
depends on the number of qubits in the chain. We found 
that for sufficiently short pulses, we could still achieve 
a fidelity of 0.99 and above for chains of up to = 9 
qubits. However, the pulse length is physically limited 
by the energy gap to higher excited states, which should 
not be excited by the pulse generator. As a consequence, 
the coupling strengths between the qubits may need to be 
reduced, which in turn increases the gate implementation 
times and could cause problems with decoherence. 

In order to achieve true scalability, error correction will 
be required. De and Pryadko recently demonstrated how 
a universal set of quantum gates could be implemented on 
a qubit lattice with Ising couplings and then implemented 
the toric code on top of this lattice to achieve scalability 
[55l |. We believe that this approach could be adopted in 
principle for our model. 
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